Helicity fluctuation, generation of linking 
number and effect on resistivity 

F. Spineanu, M. Vlad 
Association EURATOM-MEC Romania 
NILPRP MG-36, Magurele, Bucharest, Romania 
spineanu@ifin.nipne.ro 

February 2, 2008 

Abstract 

The energy of the stochastic magnetic field is bounded from below 
by a topological quantity expressing the degree of linkage of the field 
lines. When the bound is saturated one can assume that the storage 
of a certain magnetic energy requires a minimal degree of topologi- 
cal complexity. It is then possible to infer a connection between the 
helicity content and the average curvature of the magnetic field lines. 
The random curvature induce random drifts leading to an additional 
dissipation and modified resistivity. 

1 Introduction 

When the Chirikov criterion is verified for several chains of magnetic is- 
lands (developing at closely neighbor resonant magnetic surfaces in a volume 
of plasma) the magnetic field becomes stochastic. In general the magnetic 
stochasticity is taken into account in transport processes due to the high 
efficiency of energy spreading through the stochastic region. However from 
the point of view of the structure of the magnetic field it is difficult to say 
anything more than we know from Hamiltonian chaotic systems: there is 
a stochastic (exponential) instability, local Lyapunov exponents and Kol- 
mogorov length and the test particles move diffusively or have various sub- 
and supra-diffusive behaviors. 

However in a stochastic region the field must still obey some constraints. 
These constraints arise from the relation between the energy stored in the 
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magnetic field and the topological complexity of the field. The constraints 
can be briefly expressed in this way: it is not possible to support a certain 
energy in a volume spanned by transiently stochastic magnetic field lines if 
these magnetic field lines do not have a certain minimal degree of topological 
complexity. 

This should be seen in relation with the equation that expresses the topo- 
logical link in terms of writhe and twist and in relation with the dynamics of 
a twisted flux tube. If an initial amount of link is stored exclusively as twist, 
then beyond a certain level of twist the flux tube deforms and acquires writhe, 
thus distributing the higher amount of link into the two kinds of topologi- 
cal deformations: twist and writhe. In a plasma free from strong magnetic 
background (as in astrophysical plasma or solar corona) generation of writhe 
means a coiling or super-coiling instability, a large spatial deformation. In a 
tokamak the stochastic flux tubes are also subject to the writhing instabil- 
ity when a local fluctuation of the parallel electric field occurs, but they are 
more constraint by the confining B and cannot perform large spatial dis- 
placements. Instead, as a result of small deformations originating from local 
writhing (coiling) instability, they will reconnect such that, from elements 
of tubes, effectively new strings are created, with a new effective entangle- 
ment. It is reasonable to assume that these new, episodic, flux tubes, by 
their mutual linking, satisfy on the average the energy-topology constraints. 
Therefore we will assume that the field flux tubes inside the stochastic re- 
gion will reconnect to generate transiently configurations that exhibit a cer- 
tain topological entanglement. Together with the dynamical nature of the 
stochasticity phenomena, the formation of these entangled structures is tran- 
sient and we may suppose that the higher topological content results from 
a statistical average. At large time the topological reduction occurs with 
suppression of relative linking via tube merging, a process called by Parker 
topological dissipation pp. 



2 Energy and topology of divergenceless vec- 
tor fields 



For two curves 71 and 72, the link invariant is given by the formula (Gauss) 

r 

'71 J~t2 \ x ~~ y\ 



Lk ( 7l , 72 ) = ^-<f dx»<[ dy v e^ p ^-^ (1) 



This is an integer number and represents the relative entanglement of two 
magnetic lines in the stochastic region. If a line closes in itself (as a magnetic 
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line on a resonant surface in tokamak), the formula can still be applied, giving 
the self-linking. It is obtained by taking 7l = 7, 72 = 7 + en with n a versor 
perpendicular to the tangent of 7, and taking the limit e — * (this operation 
is called framing). However for a flux tube a more complex situation arises: 
the magnetic field in the tube has the lines twisted relative to the axis and 
the topological description is given in terms of the twist invariant. It is 
calculated by considering a line on the surface of the tube and the axis of 
the tube and defining the vectors: T (s), the tangent to the axis of the tube; 
U (s), the versor from the current point on the axis toward the current point 
on the line; s is the length along the axis. Then the twist is defined as 

Tw = jfds[T(s) xU(s)]-^^ (2) 

The deformation of the flux tube of axis 7 relative to the plane is measured 
by the topological number writhe, defined as 

Wr ( 7 ) = J_ I dx » I d ^ e J^vy (3) 
47r J-y J~f \x — y\ 

While the twist measures the rate at which a line twists around the axis of 
the tube, the writhe measures the rate at which the axis of the tube is twisted 
in space. The following relation exists between the three basic topological 
numbers for a flux tube 

Lk = Tw + Wr (4) 

Instead of generalizing Lk to an arbitrary but finite number of curves 
(magnetic lines) in space, it is defined an equivalent topological quantity 
(also noted Lk) which refers this time to a vector field in the volume. Instead 
of the Gauss link (a discrete set of curves) the definition will now imply a 
continuous, field-like invariant, the Chern-Simons action, which is the total 
helicity. For a divergenceless vector field £ (velocity or magnetic field) in R 3 
the helicity is defined as 

H(£)= [ cPxfacurr 1 ^ (5) 

J M 

which is the same as the integral of v • uo or A • B over volume. Consider two 
narrow, linked flux tubes of the vector field £, 71 and 7 2 . Then 

^(0 = 2Lfc( 7l , 72 ).|/Z«a; 1 |.|/Z«a; 2 | (6) 

which shows that the link of the magnetic flux tubes is a measure of the 
magnetic helicity in the volume. The total helicity in the volume is the 
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integral of the Chern-Simons form and with adequate boundary conditions 
this is an integer number, a direct consequence of the topological nature of 
the link invariant (invariance refers here to deformations of the field £ (= B) 
that do not break and reconnect the lines) 



where Ai is the magnetic potential and Fj k is the electromagnetic stress ten- 
sor. The integrand is the Chern-Simons form, or helicity density, for the 
magnetic field. The integer Q is called the Hopf invariant. The magnetic 
lines of a field B or the streamlines of a flow v are obtained form equations 
like Y^dxi/Bi = 0, generally difficult to solve. Therefore using these solu- 
tions to construct topological invariants is very difficult and we would need 
a different representation of the lines for easier handling. This is provided 
by the Skyrme-Faddeev model, or the modified O (3) nonlinear sigma model, 
where a line is a topological soliton, clearly exhibiting topological proper- 
ties (see Ward [2]). It is very suggestive that this model has recently been 
derived precisely starting from the plasma of electrons and ions, coupled to 
electromagnetic field (Faddeev and Niemi [3]). It is then legitimate to use the 
general results derived for the Skyrme-Faddeev model and in particular the 
following lower bound for the energy of the magnetic field. The inequality is 



where rj is a constant. It means that the energy is bounded from below by 
the 3/4-th power of the total helicity content in the volume or by a quantity 
that contains the total linking of magnetic lines in the volume, at the power 



A more practical measure of the topological content is the average crossing 
number C, obtained for a pair of lines by summing the signed intersections 
in the plane-projection of the spatial curves, averaged over all directions of 
projection. It differs of Eq.© by taking the absolute value of the mixed 
product. Friedmann and He [I] have extended the concept for a continuous 
field. We follow the argument of Berger [5] to find the energy bound E > 
const C 2 . 

Consider magnetic flux tubes whose ends are tied to points situated in 
two parallel planes (at distance L) and are linked one with the others. Taking 
two points on their axis, we connect them with a line and measure the angle 
formed by this line with a fixed direction in one of the plane of projection. 
This quantity is 612 . The crossing number can be expressed using this angle 




(7) 



E > 7]Q 3/4 
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3/4. 




(9) 
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The magnetic fields has magnitudes B zl and B z2 . We combine energetic 
and topological quantities by weighting the angle variation along z with the 
two magnetic fluxes 



1 



L 







C — — I dz I d x 1 I d x 2 B zl B z2 



d9 
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dz 



(10) 



The angle #i 2 is generated by the deviation of the lines (with tangent 
versors 111,2) relative to a reference straight vertical line normal to the end 
planes. This deviation is produced by the component of the magnetic field 
which is perpendicular on the main field B Z) dx.i/dz = B ±1 (x^ / B . Then 
one finds 

dOu 1 ,^ ^ N ^ , , 

-j- = — (n 2 - ni) • e 012 (11) 
dz r 12 

The versor of the line connecting points on the two magnetic field lines is 
e 9l2 =e z x r 12 . Then 

dC f f t2 t2 B Z \B Z 2 



d xid x 2 — |(n 2 - Si) • e 6l2 \ (12) 

dz J J 2nri2 

The energy of the magnetic field in a volume is Ef = ^ f d 3 xh 2 where 
b = (b x , by) = Bj_/_B 2 . By successive bounds Berger finds the inequality 

E f > const C 2 (13) 

We now have two inequalities implying the energy of the magnetic field 
in a region and two measures of the topological content in that volume: one 
is the linking of the magnetic lines, or equivalently, the helicity H, and the 
other is a more geometrical characterization of the entanglement, the average 
crossing number, C. 

Taking the inequalities as saturated we can estimate the total average 
crossing number for a certain amount of helicity H in the volume 

C ~ # 3/8 (14) 



The next step is to connect the crossing number with geometrical proper- 
ties of a generic magnetic field line. For a magnetic field line 7 the topological 
quantity crossing number C can be estimated from the number Q (7) of in- 
tersections of the line 7 with an arbitrary surface. There is a general theorem 
that allows to estimate this number as 

2irtt (7) ~ AK Y + 3K 2 (15) 
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where K i 2 are integrals of Frenet curvatures. Taking \K\ as the upper esti- 
mation of the local value of the curvature along 7, we have 

C ~ Q (7) ~ L \K\ (16) 

where L is the length of the tube. Then we have 

\K\ ~ H 3/8 (17) 

These very qualitative estimations led us to a scaling law connecting the 
average of the local curvature of a typical magnetic field line with the helicity 
inside the finite volume of the stochastic region. A magnetic line is curved 
since it is linked with other lines, and this link is generated for the magnetic 
structure to be able to store the energy in a stochastic region. Generation of 
linking also occurs when a certain amount of helicity is injected in a plasma 
volume. It is reasonable to assume that the curvature is distributed randomly 
in the volume. 



3 Effects of topology on resistivity and diffu- 
sion 

The curvature of magnetic flux tubes induces drifts of particles. Electrons 
and ions flowing along curved magnetic lines will have opposite drifts and 
local charge separations produce random transversal electric fields. For a 
finite collisionality this is a source of additional dissipation. The equation 
for ions is _ 

dfi ~ df-Q ~ 
evuEz— 1 + ev di E_ L — L - = -ufc (18) 

" OE OS 

for a transversal field E±_ ~ ^0^07!" (1 + r ) \K\ an d \K\ = |(n- V)ii|, for 
which we can use an estimation based on Eq. (1171) . We finally obtain an 
estimation of the negative current perturbation due to curvature 



8j\ 
1/2 



Von 2 T e [T~ 2 1 + p + p 2 

~e ET\ — (1 + t )\K\ — — 19 

v B l \ m,i 1 + p 



with p = j • This is not a substantial modification of the equilib- 

rium current (less than 1%), which means that the enhanced resistivity is 
mainly due to other processes. 

In general it is assumed that the magnetic reconnection does not affect 
the total helicity. However there is a dynamic redistribution of helicity (with 
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overall conservation) since in the stochastic region there are filaments of 
current and local increase of the parallel electric field, from which we have 
dh/dt = — 2E • B. Then this mechanism is a potential feedback loop: higher 
resistivity leads to higher reconnection rate and higher helicity perturbation, 
which in turn creates magnetic linking and curvature. 

In an alternative approach to the problem of topology of the magnetic 
field in a stochastic region, we can base our estimations on the much simpler 
assumption, that a magnetic line is randomly "scattered" at equal space 
(z) intervals and performs a random walk. Actually this is the classical 
assumption for the diffusion in stochastic magnetic fields. To characterize 
quantitatively the topology of the line we use the analogy with the polymer 
entanglement. A functional integral formalism can be developed (Tanaka [B]) 
taking as action the free random walk, with the constraint that the line has a 
fixed, m, winding number around a reference straight line. The mean square 
dispersion of the winding number (linking) can then be calculated 

(mVi-(lniV) 2 (20) 

where N is the number of steps. Since k\\ is low for magnetic turbulence, the 
winding number is a small number. But this represents the random winding 
naturally occuring in an unconstrained random walk of the magnetic line, 
when the magnetic perturbation is a Gaussian noise. Actually, we know that 
a given amount of helicity can only be realised by a certain volume- averaged 
mutual linking and this is an effective constraint which can only be realised 
through a much higher density of y/ (m 2 ) than the free random motion. Then 
the higher winding leads to a sort of trapping for the magnetic lines and the 
effective diffusion will be smaller that for the brownian case [7]. 

In general the use of topological quantities can improve the description of 
stochastic magnetic fields, e.g. diffusion, Kolmogorov length, etc. These are 
usually expressed in terms of mean square amplitude of the perturbation, 



but including the topological quantities can lead to more refined 



models. 

Consider the line 7 of a perturbed magnetic field and the equation Du = 
where D is the covariant derivative, similar to the velocity, v = p — A which is 
applied on a function u along the line 7. The equation says that the covariant 
derivative along the line, of the function u is zero. Then - iA) u = 

leads to it = uq exp (i <f^ dsA (s)\ . It is natural to make a generalization of 
the two-dimensional concept of point-like vortices and introduce the spinors 
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along the magnetic line 7. By the same arguments (Spineanu and Vlad [S]) 
we will need the dual (dotted-indices) spinors and we need to represent A in 
SU (2). Then, more generally, the expression of u is 



W 7 = Tr R Pexp ( i (p dx^A^ ) (21) 

7 



the trace is over the representation R of SU (2) and P is the ordering along 
7. Being a closed path this number is a functional of 7 and of A. This is the 
Wilson loop. 

We subject the fluctuations of the potential A to the constraint of min- 
imum helicity in the volume, because lower helicity allows lower energy ac- 
cording to the bound Eq.fjHJ). The Boltzmann weight in the partition function 
is then the exponential of an action representing the total helicity, i.e. the 
integral of the Chern-Simons density (compared with Eq. (17]) , here A is a 
matrix) 

S =^j M3 d 3 rs^ p Ti (a^A, - \A^A v A)j (22) 
The average of W 7 is 

j D [A] exp (S) (23) 

can be expanded in powers n~ n . The first significant term (order k" 1 ) is the 
integral of the two-point correlation of the fluctuating potential 



Tr 



{R a R b ) j dx» J X dy v (A a u (y)4(a:)) 



Wr (7) (24) 



where Wr (7) is the writhe number of the curve 7 (a, b are labels in SU (2)). 
Calculating Lfc=self-linking of the curve 7, the classical relation is obtained 
between the link, the twist and the writhe Lk (7) = Tw (7) + Wr (7). 

Therefore if the fluctuation of the poloidal flux function ip, or the z- 
component of the magnetic potential A z are such that higher helicity states 
are difficult to access, then the two-point correlations of the perturbed po- 
tential along a curve 7 can be expressed as the kernel of the Gauss integral 
for the self-linking number of 7. The result 

{ip {x\) ifi {X2)) ~ integrand of Lk (7) (25) 

possibly sheds a new light on the correlations of fluctuating magnetic quan- 
tities, since we now express them also by topological quantities. 
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4 Discussion 



We have examined the topological constraints on the stochastic magnetic 
configuration when a transient increase of helicity occurs in a finite plasma 
volume. Via bounds related to the magnetic energy that can be safely stored 
in that volume (i.e. a statistical stationarity can be attained) a scaling 
can be derived between the helicity and the average curvature of a generic 
magnetic line in the volume. The particles' curvature-drift-induced new dis- 
sipation appears to not modify substantially the resistivity. However the 
new instruments that imply the topology of magnetic field are useful : the 
average dispersion of the winding of a line relative to a reference axis serves 
to quantify the trapping of a line and the reduction of the classical magnetic 
diffusion. 

We note finally that using these analytical instruments the topological 
dissipation process may be described by the coupling of the magnetic helicity 
density (the Chern-Simons Lagrangian density) with a pseudoscalar field. 
The dynamics of this field is that of the kinematic helicity of the plasma and 
again a field-theoretical description appears to be possible. 
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